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Abstract 

A theoretical framework for dislocation dynamics in quasicrystals is pro- 
vided according to the continuum theory of dislocations. Firstly, we present 
the fundamental theory for moving dislocations in quasicrystals giving the 
dislocation density tensors and introducing the dislocation current ten- 
sors for the phonon and phason fields, including the Bianchi identities. 
Next, we give the equations of motion for the incompatible elastodynam- 
ics as well as for the incompatible elasto-hydrodynamics of quasicrystals. 
We continue with the derivation of the balance law of pseudomomentum 
thereby obtaining the generalized forms of the Eshelby stress tensor, the 
pseudomomentum vector, the dynamical Peach-Koehler force density and 
the Cherepanov force density for quasicrystals. The form of the dynamical 
Peach-Koehler force for a straight dislocation is obtained as well. More- 
over, we deduce the balance law of energy that gives rise to the generalized 
forms of the field intensity vector and the elastic power density of qua- 
sicrystals. The above balance laws are produced for both models. The 
differences between the two models and their consequences are revealed. 
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The influences of the phason fields as well as of the dynamical terms are 
also discussed. 



1 Introduction 



Since the discovery of quasicrystals by Shechtman et oOlL specimens of macro- 
scopic size have been subjected to plastic deformation [2|, y|, |4j. Dislocation lines 
have been observed by means of the transmission electron microscope j^l, 0], the 
microstructure of deformed quasicrystals being similar to the microstructure of 
deformed crystals. The first direct evidence for dislocation motion in quasicrystals 
was observed by an in-situ straining experiment by means of electron microscope 
during the study of plastic deformation of icosahedral Al-Pd-Mn single quasicrys- 
tals by Wollgarten et al . Any theory of work hardening or work softening needs 
an understanding of the interaction between dislocations, and of dislocations with 
the applied stress. For crystals the essential ingredient is the force caused by a 
dislocation in the presence of stress and it is the well-known Peach-Koehler force 
jsf . As known in configurational mechanics (jl the Peach-Koehler force density is 
the source term of the Eshelby stress tensor [10] . Obviously, a general expression 
for this configurational force caused by dislocations in quasicrystals is needed. 

In the literature, there exist different versions of generalized linear elasticity 
theory of quasicrystals. The difference of these versions lies in the dynamics of 
phonon and phason fields. Bak Hi, [12] argued that phason modes represent 



structural disorder or structural fluctuations. Following Bak's arguments in or- 
der to describe the dynamics of phonons and phasons, Ding et al 13| and Hu 
et al [lH implied that both phonons and phasons represent wave propagations. 
For that reason, the equations of motion given by them are of wave-type for 
both, phonons and phasons, assuming furthermore that the phonons and pha- 
sons have the same mass density. Such a model is often called elastodynamics of 
quasicrystals. Applications of elastodynamics of quasicrystals are given by Fan 



15l | . Fan and Mai [161 ]. and Li and Liu [17] for the study of elastodynamics, spe- 
cific heat and thermodynamic functions. Rochal and Lorman 18|, ll9] generalized 
the elastodynamics of quasicrystals in that way that the phason should possess 
a generalized phason density instead of the usual mass density of phonons. On 

2(1 , the dynamics of phasons repre- 



20] is usually called hydrodynamics 



the other hand, according to Lubensky et al 
sents a type of diffusion. Lubensky's model 
of quasicrystals. However, according to this model, the dynamics of the phonon 
fields is also of diffusive character, possessing in this way at least the drawback 
that classical elastodynamics cannot be recovered as a limit. A combination 
of the elastodynamics and the hydrodynamics of quasicrystals leads to a model 
which is called minimal model of the phonon-phason elastodynamics by Rochal 



and Lorman [19] and elasto-hydrodynamic model of quasicrystals by Fan et al [21 



According to this model, that seems more appropriate for the description of the 
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behavior of quasicrystals, the equations of motion for the phonons are equations 
of wave-type and for the phasons are equations of diffusion- type. In addition, 
several experimental results support this model as pointed out by Rochal and 



Lorman 19| . Using coherent x-ray scattering, Francoual et al [22[ have presented 
the first measurement of collective phason dynamics in quasicrystals, demon- 
strating that phason fluctuations are collective diffusive excitations. Zhu and 
Fan 23j used the elasto-hydrodynamic model for dynamic crack propagation in 
two-dimensional decagonal quasicrystals. However, as it is mentioned by Rochal 
and Lorman [3] , there exists a class of incommensurate structures where the col- 
lective phason modes are similar to sound ones and the elastodynamic model can 
be applied. Thus, both the elastodynamic and the elasto-hydrodynamic model 
have physically realistic applications. For that reason, we will deal with the elas- 
todynamic as well as the elasto-hydrodynamic model of quasicrystals in the present 
paper. 

A review, mainly from the experimental point of view, concerning electronic, 
transport, magnetic and mechanical properties of icosahedral and decagonal qua- 
sicrystals has been given by Takeuchi [24J. For the last twenty years, the elas- 
ticity theory of dislocations in quasicrystals has been an attractive research field 



25j, |26|, |27|, |28|, l29( . For an overview on the (static) elasticity theory of dislocations 



in quasicrystals we refer the reader to the articles of Hu et al 141 ] . Edagawa 30 



and Edagawa and Takeuchi [3J| . The temperature is of fundamental importance 



in the mechanical behavior of quasicrystals as it was pointed out by Bresson 



2J. Edagawa and Takeuchi [3JJ report that generally quasicrystals are brittle 
at room temperature and dislocations within them are almost immobile. They 
can only be plastically deformed at high temperatures above about 0.8 T m (T m : 
melting temperature). At sufficiently high temperature when atomic mobility by 
diffusion is faster than the dislocation velocity, perfect dislocations can migrate 
accompanying both phonon and phason strains. 

The general expressions for the displacement fields (static case) induced by 
dislocations in quasicrystals have been calculated by Ding et al 32]. Also, Qin 
et al [33j have given the analytical expressions of the displacement fields (static 
case) induced by straight dislocations in decagonal, octagonal and dodecagonal 
quasicrystals. As far as dynamics of dislocations in quasicrystals is concerned, 
using the elastodynamic model, Fan et al [34j have found the analytical expres- 
sions for displacement and stress fields induced by a moving screw dislocation in 
a one- dimensional hexagonal quasicrystal as well as the form of the energy. Re- 
cently, using the elasto-hydrodynamic model of quasicrystals, Fan et al [2l| have 
given analytical solutions for a moving screw dislocation in some quasicrystalline 
systems. It should be noted that in the mentioned works the plastic fields, e.g., 
plastic distortion, have been neglected in the equations of motion for dislocations 
in quasicrystals. However, plastic deformation is possible through dislocation 
motion at high temperatures, as observed by Rosenfeld et al [35] and Bartsch et 
al 36j|. Moreover, Li and Fan [37| have derived the Peach-Koehler force for a 
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straight dislocation in the static case for quasicrystals. In the present paper, we 
derive the general expression of the Peach-Koehler force for a moving dislocation 
in quasicrystals for both models. 



Based on compatible elasticity theory of quasicrystals, Fan and Mai [38j have 
defined the generalized Eshelby stress tensor for quasicrystals. On the other hand, 
conservation laws in physical and material spaces for a decagonal quasicrystal 



in compatible elastodynamics have been derived by Shi [39[. Shi [40[ has also 
derived the balance laws of pseudomomentum and scalar moment of momentum 
for an inhomogeneous material in the absence of body forces in the framework 
of the compatible elastodynamics of quasicrystals. Here, we find the balance 
laws of pseudomomentum and energy in incompatible elastodynamics as well as 
in incompatible elasto-hydrodynamics in presence of a continuous distribution of 
dislocations in quasicrystals. 

The paper is organized as follows: in Section 2 we set up the basic framework 
of moving dislocations in quasicrystals. Particularly, in the first subsection we 
give the geometric quantities of elastoplasticity of quasicrystals. We start with 
the compatibility conditions for the total distortion tensors and the total veloc- 
ity vectors. The terms breaking the compatibility conditions lead us naturally 
to the definitions of the dislocation density and the dislocation current tensors 
in terms of the elastic and plastic fields, fulfilling the Bianchi identities. In the 
next subsection, we give the equations of motion for the elastodynamic as well 
as the elasto-hydrodynamic model of dislocations in quasicrystals. In the third 
section, we derive the balance laws that correspond to the infinitesimal varia- 
tions of space and time, that is the balances of pseudomomentum and energy, 
respectively. The derivation of the dynamical Peach-Koehler force for quasicrys- 
tals possesses a prominent role. The above balance laws are obtained for both 
models. Eventually, the results are discussed in the last section and they are put 
into a broader context. 



2 The basic framework 

2.1 Incompatible elasticity theory of quasicrystals 

An (n— c?)-dimensional quasicrystal is defined as the projection of an n-dimensional 
periodic structure to (i-dimensional space (n > d). The n- dimensional space 
can be decomposed into the direct sum of two (orthogonal) subspaces, the d- 
dimensional physical or parallel space E\\ and the (n — aQ-dimensional perpen- 
dicular space E±, that is, K n = En © E±. In the generalized elasticity theory of 
quasicrystals there are two types of fields, the phonon (conventional) field and the 
phason field. Quantities which are associated with En are called phonon- and will 
be denoted by (■)" and quantities associated with E± are called phason- and will 
be denoted by (•) ± . Apparently, the definitions of the physical quantities in qua- 
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sicrystals are an extension of the definitions of the physical quantities in a usual 
medium enriched by the phason fields. It is important to say that all quantities 
depend only on the physical space coordinates x G E\\ . Thus, in addition to the 
usual phonon displacement field uj (x, t) G Eh , there exists the phason displace- 
ment field "U^(x, t) G Ej_. Phason displacement leads to local rearrangements of 



atoms, which are called phason flips and have a diffusive character [31| ■ 

The displacement vector U in the hyperspace E\\ © E± is given in terms of 
the phonon and phason displacement fields as follows 

U=(u\u ± ) G E\ } ®E ± . (1) 

The total distortion tensors, which are the gradients of the components of the 
displacement U with respect to the spatial physical coordinates, and that corre- 
spond to the fields of phonon, M- T , and phason, /3^ T , can be decomposed into 

elastic /?]'., /3^j and plastic parts , f3£j F ', respectively 

ul :=4 T = 4 + 4 P eE\\®E h (2) 
<j ■= Af = fij + ^i® #|] , (3) 

where comma denotes differentiation with respect to the spatial physical coor- 
dinates. Similarly, we can decompose the time derivative of the displacement 
fields 

uj := v } T = v }+v} P EE h (4) 
u^=vr = v^ + vr eE ± , (5) 

where v\ and are the elastic velocities corresponding to the elastic fields of 

IIP IP 

phonon and phason, respectively and v- and vf are the plastic or initial ve- 



locities [41[ in the spaces E\\ and E±, respectively. A superimposed dot denotes 
differentiation with respect to time. 

Due to the definition of the displacement field U in the hyperspace (see 
eq. ([I])), the compatibility as well as the incompatibility conditions will be analo- 
gous for the phonon and phason quantities. Therefore, the total fields satisfy the 
compatibility conditions 

e jkiP!i T k = 0, ZjkiPuJ = 0, (6) 

C - 4J = 0, Af - vff = 0, (7) 

where e^-jy is the permutation tensor or Levi-Civita tensor. On the other hand, the 
elastic as well as the plastic fields satisfy separately the following incompatibility 
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conditions 



and 



4 = e i*4,* e^i®^ih (8) 

ai = e 3M /3u, k eE ± ®E h (9) 

4 = 4-^ GEn®% (io) 

= eE ± ®S||, (li) 



= - e iH/3£, (12) 

3-LP 



= - e mPil,ki ( 13 ) 

4 = -4 p + u !J> ( 14 ) 
% = ( 15 ) 

The first two incompatibility conditions define t/je phonon and t/ie phason dislo- 
cation densities, ac-j and a^, respectively. The last two incompatibility conditions 

define the phonon and the phason dislocation current tensors, ifj and 1^ , respec- 
tively. The dislocation current tensor describes the movement of dislocations and 
contains rate terms. One can read more about the significance and the phys- 



ical interpretation of the dislocation current tensor in [42|. At this point it is 
necessary to generalize the well-known Bianchi identities 42|, |43) for the disloca- 
tion density and the dislocation current tensors to quasicrystals. The dislocation 
density and the dislocation current tensors fulfill the Bianchi identities 

°L = °> = °> ( 16 ) 

«!j = e 3ki4, k , ®ij = £jkiiu, k - (17) 

In a quasicrystal, a (perfect) dislocation is a line defect with the following 
Burgers vector B in the hyperspace En © E± 

B = (b^tr 1 ) G E\\@E ± , (18) 

where 

6| = jf4^ and ^ = jf/^ Xj (19) 

are the components of Bi in the spaces En and E±, respectively, and C is a 
closed contour surrounding the core of the dislocation in the physical space E\\. 
Alternatively, one can express the Burgers vectors in terms of the dislocation 
density tensors using the Stokes theorem 

h \= I and h t = I OySj, (20) 

Js Js 
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where S is the surface bounded by the curve C and pierced by the dislocation. If 
we integrate the Bianchi identities ( ITT]) over the surface S and we use the Stokes 
theorem and the above formulas for the Burgers vectors, we obtain the following 
formulas 



\ = f c 4* and T = f c ^^ (21) 



which are actually conservation laws for the Burgers vectors. The integrals on 
the right-hand side give the flux of the Burgers vectors through the contour C 
per unit time and for that reason the dislocation current tensor is also called 



dislocation flux tensor [42 



2.2 Equations of motion of dislocations in quasicrystals 



We start our study with the elastodynamic model of quasicrystals. The La- 
grangian density C is given in terms of the kinetic energy density T and the 
elastic energy density W as follows 



C = T — W 



with 



T 



~-A v l 



1 



;Pcfi {VI 



(22) 



(23) 



where p is the usual mass density and p e g is the effective phason density [18|, [19 



For the unlocked stated of a quasicrystal, that means for T > Tc (Tc- transition 
temperature) in 3D-systems 44|, the elastic energy density W can be written in 
a quadratic form of phonon and phason strains 



(24) 



where C^i are the elastic constants in the classical elasticity theory, are 
the elastic constants of the phason field in E± and are the elastic constants 
associated with the phonon-phason coupling with the following symmetries 



a 



ijkl 



klij 



a 



ijlk 



jikl i 



D 



ijkl 



Da 



jikl ; 



E 



ijkl 



E, 



klij • 



(25) 



The specific form of the tensors of the elastic constants and of the elastic energy 
density depends on the considered type of quasicrystals and has been examined 
for various types of quasicrystals for example pentagonal and icosahedral 46 



octagonal and dodecagonal [47j, cubic [48| as well as for two-dimensional qua- 
sicrystals [49| . 



1 For the locked state there is a discussion whether the elastic energy density can have a 
quadratic form or the phasonic elastic energy density should depend linearly on the phason 
strains 0, HH, Q Efl . 
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Here, the Lagrangian density C is considered as a smooth function of v" , v- 1 , (3" 
and f3 ± . The volume integral of C is given by 

£(vll,v ± 5/ 3ll, / 3 ± )x ) (26) 



n 

where Q is an open, connected subset of R 3 (d = 3) with smooth boundary 
dfl = S. The constitutive relations for the considered Lagrangian density are 

ii dC ii / v 

pf = Tnr = wl l , (27) 

^ = = P& V i ' ( 28 ) 

II dW II ± 

± 9iy ii ± 

a ij = = D klijPkl + E wPkl- ( 30 ) 

Evidently, p\ and p^ 1 are i/ie phonon and phason momentum vectors and cr|- and 
<T^ are the phonon and phason stress tensors, respectively. a£j describes the stress 
components along the Xi direction in E 1 - applied on the surface orthogonal to the 
Xj direction in EK Notice that the phonon stress is symmetric, afj = a^, but 
the phason stress is not, ^ aji, since the relevant constants D^i and E^i 
lack the necessary symmetries (see also, e.g., [H]). Furthermore, the kinetic and 
the elastic energy densities, using the constitutive relations (|2"7|) - (|3"U|) . can be 
expressed as quadratic form 

T = \{p\v}+pivi) ) (31) 

^ = ^(44+Ai4)- (32) 

Introducing a generalized (phason) body force density and the conventional 
(phonon) body force density /J', the equations of motion can be written in the 
form 

P i i-°'iij = fi> ( 33 ) 

Pi - ojjj = ft- (34) 

By substituting the constitutive relations (l271)-(l30l) in the above equations of 
motion and making use of the relations ([2])-([5]), we obtain the inhomogeneous 
partial differential equations for the displacement fields u" and u- 1 in a homoge- 
neous material 

puj ~ Qjkiuljj - AjW«fc,y = pv l } F ~ Ci jk iPlf d ~ DijklPuj + f}, (35) 
p cff ui - E^iui^ - Dkujul^ = p eS v^ p - D mj pl* - Ew^ + ft, (36) 
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where the plastic fields and the external forces are given as source terms. Eqs. (|35|) 
and (136|) are wave- type equations for the displacement fields u" and u -1 , respec- 
tively, for anisotropic media and they are coupled due to the constitutive tensor 
Dijki. Eqs. ( 1351) and (1361) represent the equations of motion for the incompatible 
elastodynamics of quasicrystals. In the case that = p we recover the model 



that was proposed by Ding et al [13| and Hu et al [1 



We continue our study with the elasto-hydrodynamic model of quasicrystals. 



This model 19] introduces the attenuation of phason modes which is essential in 
the dynamical theory assuming that the phason bulk forces ff- = —Du^-, where 
D is the friction coefficient. In addition, the phason momentum in quasicrystals is 
not conserved implying p eS = 0. Therefore, the phasons do not give a contribution 
to the kinetic energy. Under these circumstances, neglecting the phonon body 
force density ff = and putting for the incompatible case f^ = —Dv^, the 
equations of motions (133 p and (I34p take the form 

ojjj = Dvt. (38) 
The above equations can be written more analytically as follows 

f>4 ~ C ijkiA,ij ~ D W u k,u = p4 P - CijkiPklj ~ DijkiPuj, (39) 
Dili ~ E ijkiui :lj - D klij ul tlj = Dv^ v - DuypHj ~ E ijkiPii V j- (40) 

The first eq. ( 139]) is a wave-type one for the displacement field u" while the second 
one (eq. ( HUl) ) is a diffusion-type equation for the u^, both for anisotropic media. 
Again the plastic fields are the source terms for the phonon and phason displace- 
ments and of course the two equations are coupled. Eqs. (|39|) and (f40|) are the 

equations of motion for the incompatible elasto-hydrodynamics of quasicrystals. 



3 Translational balance laws 

In this section, we derive the balance laws that correspond to the infinitesimal 
variations of space and time. It is clear that the translations in space will give 
the balance law of pseudomomentum while the variation of time will lead to the 
balance law of energy. We consider an arbitrary infinitesimal variation of the 
fields, leading to an infinitesimal variation of I 

61= [ 5£x= f {5T-5W)*. (41) 
Jn Jn 
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3.1 Balance law of pseudomomentum— Peach-Koehler force 

We specify the variation here to the translational variation in the space coordi- 
nates as 

* = ^A (42) 

where ^f— is the generator of the translation group in space and e m are the 
parameters of the translation group. Then, eq. ( 14"T|) reduces to 

f / dT dW\ 

61 / b^-^K*- (43) 

Jn X( Jx m ox m / 

We begin again our study with the elastodynamic model. Firstly, using the 
incompatibility conditions ( flOj) and (llip . the calculation of the gradient of the 
kinetic energy density gives 

Ox 9t ' Pi Pirn) Pi Pirn, Pi Pirn Pi^im Pi ^im' ( 44 ) 

Also, the gradient of W, using the constitutive relations ff29l) and fl30l) and the 
symmetries of the coefficients Cy&j and E^m, gives 

<9J¥ || || ± ± 

~pT~ = a klPkl,m + a klfikl,m- (45) 

Multiplying the compatibility conditions (jSJ) and (j5J) with the permutation tensor, 
we obtain two auxiliary identities which will be used in further calculations 

P\l,m ~ @km,l = e mlj(^kji 

Pkl,m ~ Pkm,l = e mljCXkj- (46) 

Next, if we substitute the relations (j4"6"|) in eq. ( I4"5j) . we get 
dW 



8Xr 



a kl e 



mlj^kj + a klPkm,l + a kl e mlj a kj + ^hAmd' ( 47 ) 



Inserting the eqs. (|44p and f|4"T|) into eq. (j4~3l) and rearranging the terms in order 
to use the equations of motion fl33l) - fl34l) . we obtain the following 



(4* 



- PUL ~ Pt J L - e mlAAj ~ e mlj^a a ij ~ fi Pirn ~ fi~PL) e m*- 

On the other hand, we can also write 



51 = ( 5£x= ( J^-e m x= / -^-(£5 m i)s m x. 
Jn Jn ox m J Q ox\ 



(49) 
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Combining eqs. (HHj) and (|49|) . we arrive at a translational balance law for qua- 
sicrystals in global form 

— (-£ 5 mZ - o-JI^L ~ °« A>J + ^(pI^L + pIPL))* ( 50 ) 



In the above balance law we can recognize in the first integral two important 
quantities, the first one that enters in the divergence is the Eshelby stress tensor 

P m i := -£ 5 m i - al/3f m - o-afiln (51) 

and the second one that enters in the time derivative is the pseudomomentum 
vector 

V m :=-v\A m -V^L- (52) 

In elasticity, the tensor P m i is also called the static energy-momentum tensor [lo[ 
and is familiar in a general field theoretical concept 50|. The integrand of the 
right-hand side integral is constituted by the sum of two forces, the dynamical 
Peach-Koehler force density 

F m K ■= P l ! ll L + Pt J L + e 'miAAj + emijO-uaij (53) 
and the Cherepanov force density 

F° := S}fiL + ft PL (54) 

which is produced by body forces in presence of elastic distortions. Both the 
Peach-Koehler force and the Cherepanov force are configurational forces and 
they are also called material forces (9)]. Their discussion, familiar from simple 
elastic media, can be taken over. For their elastic counterparts one can see the 



dynamical Peach-Koehler force in [5l|, |52| and the Cherepanov force in [53j, 154 



After the definitions (i5"T|) - (i5"%|) . the balance law of pseudomomentum ( 150|) can be 
elegantly written in global form as 

JM L^ +FS * (55) 



or in local form 

Pip , Fi-d 

+ F%. (56) 



dP-ml 9 T'm j-iPK 



dxi dt 

It is evident from the obtained results, that all physical quantities that enter 
into the balance law of pseudomomentum are strongly affected by the phason 
component of a dislocation. The balance law of pseudomomentum has been 
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derived by Shi 40|] for inhomogeneous materials in absence of body forces in the 
framework of compatible elastodynamics of quasicrystals. 

Furthermore, from the integral form of the balance law (l55l) . we can obtain 
the dynamical J -integral 



J n 



dP m i 



dV r , 

dt 



PminS 



dt 



(57) 



^PK 



n 



where n is the outward unit normal vector to the surface S. 

We proceed to find the form of the Peach-Koehler force for a straight dis- 
location. In analogy to the classical elastodynamics of dislocations in crystals 
42l . l55| . we give the dislocation current and the dislocation density tensors for a 
straight dislocation to quasicrystals as follows 



kj = e jkiV k ot u , 



l3 -jklVkWu, 

alj(x -x',y- y') = b}TjS(x - x')5(y - y'), 
a ij( x -x',y- y') = b t T j$(x - x')8(y - y'), 



(58) 
(59) 
(60) 



where x' = (x',y f ) is the position of the dislocation line and V = x' is the 
velocity of the moving dislocation, r = r(x) G E\\ is the tangent vector to the 
dislocation line and <$(■) is the Dirac delta-function. Using the relations ([58]) the 
Peach-Koehler force density ( |53l) becomes 



-iPK 



UJ a ij 



(61) 



Further, making use of the expressions ( 1591) and 1601) . eq. ( IBTj) gives the form of 
the dynamical Peach-Koehler force for a straight dislocation in quasicrystals 



PK 




,nij (pjVi + (rl)(x', y')a.l(x -x',y- y') 
+ (Pi~V t + a^)(x', y')a±(x -x',y- y') 




(pjVt + al)(x', y')b\ + (piV t + ai)(at, y')bi r, 5(x - x')5{y - y')x'y 



(62) 



It is obvious that the generalized dynamical Peach-Koehler force of a moving 
dislocation with velocity V and Burgers vector B in quasicrystals is produced by 
the phonon and phason stresses and the phonon and phason momenta. If we omit 
the dynamical terms in eq. ( 162|) . then the obtained force is in accordance with 
the Peach-Koehler force that has been derived by Li and Fan 37J. In addition, 
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eq. fl62l) reduces to the classical dynamical Peach-Koehler force for crystals (see, 
e.g., [43|) when the phason fields are absent. 

In the case of the elasto-hydrodynamic model of dislocations, the general 
formulas of the balance law of pseudomomentum ( )55l) or (1561) are still valid. 
However, the quantities that constitute this balance law are modified as follows: 

- The Eshelby stress tensor is given by the same formula (15 ip 

P ml := -C 8 ml - <t|/C - a^PL (63) 

with a difference in the explicit form of the Lagrangian density C = T — W, 
because the kinetic energy density consists in the current case only of the 
phonon part, T = 

- The pseudomomentum vector is given by 

V m ■= -p!pL (64) 

It is obvious that only the phonon fields give a contribution to the pseudo- 
momentum vector, which has the same form as in incompatible elastody- 
namics of crystals. 

- The dynamical Peach-Koehler force density takes the following form 

F m K ■= pUL + e mijO-lal + emijO-^aiy (65) 



- The Cherepanov force density is given by 

Fm ■= -Dv^L- (66) 

The physical meaning of the Cherepanov force density is now a configu- 
rational force density caused by the phason friction force f -1 = -Dv 1 in 
presence of the elastic phason distortion (3^. Thus, the term Dv 1 influ- 
ences the J-integral only through the force density (16^1) . 

Furthermore, in the case of a straight dislocation the dynamical Peach-Koehler 
force density is given by the same formula as eq. (1621) for pf = 0. 



3.2 Balance law of energy 

We continue in this fashion to find the balance law that corresponds to the time 
translation 

f =4- (67) 
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where & is the generator of the translation group in time and e is the parameter 
of the corresponding group. The corresponding variation of X will be 

We examine first the elastodynamic model. Due to the quadratic forms of T 
and W (eqs. ( 1311 and (1321). we have respectively 

dT II II i i 

— =p» v » +p . v . (69) 

and 

dW II • II 

= 44 + 44- (70) 

Using the compatibility conditions (flOj) and (llip . eq. (!70|l can be further written 
dW 



dt 



Inserting eqs. ( R)9"j) and ( 17 ip in the expression ( l68p and rearranging the terms in 
order to use also the equations of motion (!33j) - (l34l) . we finally obtain 

61 = f (J^W + Pi v t) ~ ~^( a ij v i + 4^) - 44 ~ 44 - fi v i ~ fi ±v f) £ *- 

(72) 

On the other hand, 

61= [ 5£x= [ ^ex. (73) 
Jn ■ Jndt ■ 

Combining eqs. ( T7Z|) and (|75|) . we obtain t/ie balance of energy for quasicrystals 
in global form 

^(pW - £ ) - ^r(4 u l + 4^ 
(44 + 44 + /N 1 + z^)*- ( 74 ) 

n 

In the first integral, the quantity in the time derivative is the Hamiltonian density 

H = p\vl+pivi-C = T + W (75) 
and the quantity in the divergence is the field intensity vector 

Sj = 4 V ! + a ii v i ( 76 ) 
14 



which describes the material energy flux. Finally, in the integral on the right-hand 
side we can see the elastic power density 



A = °IA + °¥ti + fS + fi"i, ( 77 ) 

which consists of two parts, the first one describes the energy exchange between 
the moving dislocation and the quasicrystal and the second one is the power pro- 
duced by the body forces. It can be seen that the elastic power density for a 
quasicrystal is produced by phonon and phason stresses and phonon and phason 
body forces. If we compare eq. (1771) with the mechanical power in generalized 



continua (see, e.g., 56]), we observe that the stress tensors have the same dimen- 
sion as the dipolar forces and the dislocation current tensors play the role of the 
generalized velocities (velocity gradients). 

Taking into consideration the definitions (1751) -(177j). the balance law of energy 
(I74j) can be rewritten as 



V dt dx 



Ax (7S 



or in local form 

z :;; 

In the same manner as in the previous subsection, we proceed to find the 
balance law of energy according to the elasto-hydrodynamic model of dislocations. 
The general formulas of the balance law of energy (1751) or (!79|) are still hold. 
However, the quantities that take part in this law are modified as follows: 

- The Hamiltonian density changes due to the change of the kinetic energy 
density T 

U = T + W 1 T= l -p\v\. (80) 

- The field intensity vector does not change 

S^afA+vH- (8 1 ) 

- The elastic power density is modified as follows 

k = a\ 3 ll + ^-Divif. (82) 

The term Dv 1 enters the balance law ( !78l) through the last term of the 
elastic power density ( 1821) . 
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4 Conclusion 



We have embedded the dynamics of dislocations in quasicrystals in a general dy- 
namic field theoretical framework. For the quasicrystals the instruments of field 
theory must be applied to an enlarged space, because the phason space is added 
to the phonon space. Formally such space extension is akin to adding piezo- 



electricity, piezomagnetics and magnetoelectricity to elasticity (see, e.g., [571]). 
Procedurally the necessary extensions are straightforward. 

The dynamics of quasicrystals is a complex phenomenon and there is not yet 
an agreement, from the theoretical as well as from the experimental point of view, 
for the form of the generalized elasticity theory in dynamics of quasicrystals. In 
the present work, we cover two models that are mostly used in the literature 
for the description of the dynamic behavior of quasicrystals. The differences 
between these two models are essential since the governing equations for the 
elastodynamic model are of wave-type for both, phonons and phasons, while for 
the elasto-hydrodynamic model are of wave-type for the phonons and of diffusion- 
type for the phasons. The last years, the elasto-hydrodynamic model has gained 
the interest of the researchers since it seems that physically describes better the 
motion of quasicrystals. However, we should note that the equations of motion, 
corresponding to the last model, consist a coupled system of partial differential 
equations of wave-type and of diffusion-type thereby presenting some principal 
difficulties to be solved analytically. 

Furthermore, we have defined the dislocation density tensor and the disloca- 
tion current tensor for quasicrystals in terms of the phonon and phason fields. 
Using the field theoretical framework, we have introduced the Eshelby stress ten- 
sor, the pseudomomentum vector, the Hamiltonian density and the field intensity 
vector for a quasicrystal with incompatible fields. By means of these quantities 
the dynamic translational balance laws are established. As source terms of the 
translational balance laws we have obtained configurational forces, namely the 
Peach-Koehler force and the Cherepanov force in quasicrystals and the elastic 
power caused by phonon and phason stresses and phonon and phason body forces. 

Finally, we would like to remind the reader that our results are valid at suf- 
ficiently high temperatures since we deal with movement of dislocations. For 
instance, one can see the exact range of temperature that can be used in ex- 
periments as well as results from in situ straining experiments concerning the 
dynamic behavior of dislocations for icosahedral Al-Pd-Mn single quasicrystals 
in [581 . 
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